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Abstract 

We derive equations of motion for the tachyon field living on an unstable non-BPS 
D-brane in the level truncated open cubic superstring field theory in the first non-trivial 
approximation. We construct a special time dependent solution to this equation which 
describes the rolling tachyon. It starts from the perturbative vacuum and approaches 
one of stable vacua in infinite time. We investigate conserved energy functional and 
show that its different parts dominate in different stages of the evolution. We show 
that the pressure for this solution has its minimum at zero time and goes to minus 
energy at infinite time. 
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1 Introduction 

The physical problem which we would like to address in a series of papers is a time dependence 
of the tension of an unstable non-BPS D-brane in the process of its decay. It is expected 
that starting from its maximal value, equal to the density of the vacuum energy, the tension 
of an unstable non-BPS D-brane becomes zero when the tachyon field approaches the true 
vacuum. In this paper we start with study of classical time dependent solutions describing 
a motion of the tachyon on an unstable non-BPS D-brane within level truncation scheme in 
cubic open superstring field theory. In the next paper we will take into account higher level 
fields including vector fields and study a dynamics of a D-brane tension. 

The tachyon potential and various classical solutions in SFT on an unstable D-brane 
system have been intensively studied in the last two years (see [U El El I3J E] , reviews [HI E] 
and references therein) . In the beginning only time independent solutions representing either 
the tachyon vacuum or static D-branes of lower dimension have been examined. The process 
of production and decay of unstable branes can be described as a space-like D-brane jUj • For 
earlier study of the tachyon dynamics on brane-antibrane system see, for example, [TO] . 

Recently, Sen has proposed to study the classical decay process of an unstable D-brane 
in the framework of SFT [TT| IT2*| IT3*1 IT4] . As compare to the usual local field theory the time 
dependent equation of SFT containing a tachyon has special features. First, this equation 
involves infinite number of other fields. Second, the string field theory action contains 
infinite number of time derivatives, and hence the initial value problem seems to be not 
well-defined. Nevertheless, Sen [TTJ has shown that it is possible to construct a family of 
classical solutions to the string field theory equations of motion characterized by the initial 
position and velocity of the tachyon field. 
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There are at least two different ways to construct time dependent solutions within SFT: 



• One can study solutions to SFT equation of motion 



QA + A * A = 



(1.1) 



perturbatively, specifying some initial data. This way has been followed in the back- 
ground approach to scattering theory in SFT in [TS] as well as in Sen's paper [Tlj . 

• One can also study a tachyon dynamics within the level truncation scheme. As it has 
been noted in there are many similarities between the tachyon equation of motion 
obtained in SFT via the level truncation and the p-adic string equation of motion 



The goal of this paper is to study the tachyon dynamics in SSFT on a non-BPS D-brane 
within the level truncation scheme [213] . We derive the equation of motion for the tachyon in 
the first nontrivial approximation. We construct a special time dependent solution to this 
equation which starts from the perturbative vacuum and approaches one of stable vacua. To 
construct this solution we consider it as a positive time part of a solution that interpolates 
between two non-perturbative vacua. As in the case of the p-adic string [Hi] to perform 
numerical study of the SSFT tachyon equation of motion with boundary conditions on ± 
infinities, it is useful to deal with the integral form of the equation of motion. The integral 
form of the tachyon equation of motion in SSFT is similar to the integral form of the tachyon 
equation of motion in the p-adic string ^H] for p = 3 while the integral form of the bosonic 
SFT is similar to one of the p-adic string for p = 2 ^HIEE]- In contrast to p-adic case, kernels 
of the integral equations describing SFT and SSFT tachyons are not positively defined and 
application of an iteration procedure is more subtle. By our request these simulations have 
been performed in [2H] and we will use these results in the paper. Numerical analysis [2H] 
is based on an iteration procedure of solving nonlinear integral equations. It is interesting 
to note that SFT and SSFT tachyons equations can be considered as perturbations of the 
corresponding p-adic string equations of motion. The superstring interpolating tachyon 
solution can be obtained in a linear approximation around the corresponding p-adic (p=3) 
solution. This approximation rather well reproduces results of numerical calculations. 

In order to understand qualitative properties of solutions we write down the conserved 
energy functional. The energy E is the sum of three terms, namely the kinetic term Ek, the 
potential term E p and the non-local term E n i. There are no difficulties to study numerically 
the kinetic E^ and potential E p terms. To investigate numerically contribution of the non- 
local term E n i we write it in a special form admitting an integral representation. 

The pressure is also the sum of three terms 



The non-local term E n \ v is just a part of the non-local term E n \ contributing to the energy E. 
We show that for large time P(oo) = —E while for t = 0, under a special approximation for 
the tachyon solution to the SSFT equations of motion (see Sect. 2.2), -P(O) ~ E. The energy 
E is equal to the tachyon potential minimum and according to Sen's conjecture the tension 



[m MUM- 



Pit) 



E + 2E k (t)+E nlp (t). 



(1.2) 
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of the non-BPS D-brane compensates this energy and the total energy of a system is equal 
to zero. This compensation has been checked for cubic SSFT with rather good accuracy 
[2*U] . Since the total energy is equal to zero the total pressure also vanishes at large time for 
our solution. We also show that the pressure for this solution has its minimum at zero time 
and goes to minus energy at infinite time. 

Let us also note that there are other approaches to study rolling tachyons. One of them 
is related to CFT description of a non-BPS D-brane dynamics. Sen has noted [TT| IT2*1 HI] 
(and refs. therein) that dynamics of a tachyon on an unstable D-brane can be described 
by an effective Born-Infeld action with an exponentially decreasing potential. Though this 
action cannot be derived from first principles it has some features desirable in a classical 
dynamics of open SFT around the tachyon vacuum. Namely, in this effective field theory 
there are solutions with fixed energy density and asymptotically vanishing pressure. These 
solutions are also interesting from cosmological point of view, see [2J for details. 

Let us note that it would be interesting to see more relations between SFT and the 
effective Born-Infeld action. An incorporation of vector field could clarify this problem. 
This hope is supported by the observation that an incorporation of vector fields, apparently, 
reproduces the Born-Infeld action |2*2~] 123| 12*1]. 

The paper is organized as follows. In Section 2 we introduce a general setup, namely the 
cubic superstring field theory action for a description of an open superstring living on a non- 
BPS D-brane. Within the level truncation scheme we get from this action tachyon equations 
of motion. In fact in this paper we restrict ourself to the first non-trivial approximation 
which includes the tachyon field and one auxiliary field. Under an assumption that even 
derivatives of the auxiliary field vary slowly we get an approximate version of the tachyon 
equation of motion. Then we rewrite differential equations in the integral form. In Section 3 
we present results of numeric study of the equations obtained in Section 2. 

In Section 4 functionals of the energy and pressure are studied. Conservation of energy 
is calculated explicitly. For the pressure it is shown that it has minimum at zero time and 
equal to minus energy at infinite time. 



2 Equations of Motion for Tachyon Field on Non-BPS 
D-brane 



To describe the open string states living on a single non-BPS D-brane one has to consider 
GSO± states PQ. GSO— states are Grassmann even, while GSO+ states are Grassmann 
odd. The unique (up to rescaling of the fields) gauge invariant cubic action unifying GSO+ 
and GSO— sectors is 



S[A+,JL 



1 

^L2 



1 «y_ a | Ah Q b A + )) + ^«Y- 2 |Ah Ah A+)) 



+-((Y_ 2 \A_,Q B A_ 



((Y_ 2 \A + ,A_,A_)) 



(2.3) 



For A_ = one gets the action [23 123 • Here the factors before the odd brackets are fixed 
by the constraint of gauge invariance, that is specified below, and the reality of the string 
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fields A±. Variation of this action with respect to A + and A_ yields the following equations 
of motion (see [2U| for notations and details) 

QbA± + A + * A± — A_ * A T = 0. (2.4) 

The action (J2.3j) is invariant under the gauge transformations 

6A± = Q B A± + [A ± , A+], +{A T , A_}, 

where [,] ({ , }) denotes *-(anti) commutator and A± are gauge parameters. 

To perform actual calculations the level truncation scheme is employed. 

Let us denote by Hi the subset of vertex operators of ghost number 1 and picture 0, 
created by the matter stress tensor Tb, matter supercurrent Tp and the ghost fields b, c, d£, 
r\ and <fi. We restrict the string fields A+ and A- to be in this subspace 7i\. 

Next we expand A± in a basis of Lq eigenstates, and write the action ()2.3|) in terms of 
space-time component fields. The string field is now a series with each term being a vertex 
operator from 7i\ multiplied by a space-time component field. We define the level K of 
string field's component Ai to be h + 1 , where h is the conformal dimension of the vertex 
operator multiplied by Ai, i.e. by this convention the tachyon is taken to have level 1/2. To 
compose the action truncated at level (K, L) we select all the quadratic and cubic terms of 
total level not more than L for the space-time fields of levels not more than K. Since our 
action is cubic, IK < L < 3K. 

2.1 Action and equations of motion 

The level (1/2,1) action contains the tachyon field <fi and one auxiliary field u 



9o a 2 



a' 1 1 

u2 ( x ) ~ —d^(p(x)d' l (p(x) + -4> 2 {x) - —u(x)(f)(x)(j>(x) 



(2.5) 



where the tilde is the action of the operator 

Du = exp(— a' In 7 □), <f> — Du4>, etc. (2.6) 

7 = ^7^, □ = and the signature is ( — h + + . . .). The equations of motion have the 
form 

2u{x) - —D U ((D n <P) 2 ) (x) = 0, (2.7a) 

(<*'□ + 1/2) <j>{x) - ^- 2 D n ((£>□<£)(£>□«)) (x) = 0. (2.7b) 

Applying Du to both equations above we can write them as 

1 



u ~ ^D 2 n ((0) 2 j = 0, (2.8a) 
a'D + 1/2) 4> - ^D 2 D (u^j = 0. (2.8b) 
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Here D^cj) = exp(— 2a' In 7 D)0 and for simplicity we do not write arguments of functions. 
Performing a rescaling of fields and coordinates 



\[2 ~ 4 
0(7) 07^ 



with 



a 2 = 4a = —4a' In 7 and q A 



a 



We recast the latter system to the following form 

T(x) - Dq (vP 2 ) (x) = 0, 



1 



2a 4 In 7 



(2.9) 



(2.10) 



'2.11a) 



(q 2 D + 1) - £>□ (Ttf ) (2) = 0. 



Hereafter Drj = ex p(|D)- 



(2.11b) 



2.2 Approximation u = u 

The approximation it = u is valid when even derivatives of the field u vary slowly. Substi- 
tuting u = u in action ()2.5|) we are left with the action which after integration over u and 
rescaling has the form 



dx 



where 



~V> 2 (z) - q -d^{x)d^{x) - l -^\x) 



1, 



ip(x) = \l/(x) = exp(-n)i/j(x). 

8 

Equation of motion for this action is 

(q 2 D + l)^(x) = D 2 n (* 3 ) (x), 
which can be presented as a system of two equations 

T(x) - V 2 (x) = 0, 



(2.12) 



(2.13) 



(2.14) 



(2.15a) 



(q 2 D + 1) V(x) - D 2 U (Ttf ) (x) = 0. 



(2.15b) 



2.3 Integral equations for space homogeneous configurations 

For space homogeneous configurations one has 

T(t) - D 2 (^ 2 ) (t) = 0, (2.16a) 
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V^ + lJ*W-A 2 (T*)(t)=0. (2.16b) 
Initial and tilded fields are related as 

ip{t) = = e~^i;(t), (2.17a) 

v(t) = T(t) = e~**Pv(t). (2.17b) 

Note that the operator _D^ = e 4—7 is not well-defined on an arbitrary smooth function. 
This happens because the operator in the exponent is positively defined. However, we can 

_____ _„ _____ 

rewrite the system of equations (|2.16j) using only the operator D t = e 4 ^ 7 . The resulting 
system is 

D; 2 T = (2.18a) 
d 2 \ 

~ q2 lP + 1 ) D * 2 * = T *' ^ 2 ' 18b ^ 

As in the case of p-adic string [T^] one can rewrite this system of differential equations in 
the integral form 

CT = (2.19a) 

(-q 2 ^- 2 +l)CV = ry, (2.19b) 

where 

= Ci/4* 

and 

C a m = -== / e"— ^(t')dt'. (2.20) 
V47ra j_oo 

The initial fields are related with tilded ones by means of the integral operator C a 

v = Ci /B T = Cx/nv, (2.21a) 

^ = Ci/8* = C 1/8 f (2.21b) 
In the case q 2 = equations ()2.18j) take a simple form 

CT = * 2 , (2.22a) 

Cm = Ttf. (2.22b) 
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2.4 Approximate integral equation 

On the space homogeneous configurations instead of (j2.14|) one gets 



(2.23) 



This equation can be rewritten in the integral form 



# 3 , 



(2.24) 



where = Ci/ 4 \1/ and C1/4 is defined in ([2.20|h If we drop the term with the second 
derivative we are left with 



This is a p-adic effective field equation for p=3 [T71 EH] . For the bosonic string theory one 
gets equation (|2.25|) with p=2. As it has been shown in [TH] properties of p-adic equation 
are drastically different for p=2 and p=3. In particular, for p=2 there are no rolling tachyon 
solutions. 

3 Solutions to Integral Equations 
3.1 Interpolating solutions 

Equation (J2.25|) has been studied numerically in [Fi]\ for p-adic effective field equations for 
arbitrary p and recently in j!4j . More general equation ()2.24|) as well as a system of equations 
()2.19|) has been studied in [22] • Here we review results of numerical calculations performed 
in [21]. It was shown that: 

• Solution of (|2.25|) interpolating between two vacua is given by a monotonic function 
^0 presented in FigJTk. 

• Solution of (|2.24|) for q 2 < q 2 r = 1.38 is given by a function \1/ interpolating between 
two vacua (see FigJTjD). This function oscillates near with exponentially decreasing 



• Function \1/ solving the system ()2.19|) for q 2 < q 2 cr = 2.2 is an odd function with a jump 
at t = and ^(±00) = =Fl. T is even function. In Fig|2t these functions for q 2 ~ 0.96 
are presented. 

• The smoothed functions v and if), defined by (|2.21|) . are continuous, see Fig|2Jo. 

• For q 2 ~ 0.96 function ip ()2.21b|) does not differ much from the function obtained from 
the solution of (j2.24J) . by the same formula (|2.21b|) ip app = Ci/s^ app (see FiglSJ). 



Cm = W : p = 3. 



(2.25) 



amplitude. 
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Figure V. a) ^ = ^a PP for q 2 = 0; b) ^ app for g 2 = 0.96. 
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Figure 2: a) \1> (thin) and T (thick) for q 2 = 0.96; 

b) smoothed functions: ip (thin) and v (thick) for q 2 = 0.96. 




Figure 3: Smoothed functions ip app (thin) and ip (thick) of solutions for q 2 = 0.96. 
3.2 Rolling solutions 

The iterative procedure has been formulated so that the function ^ is an odd one and the 
function T is even. Results of calculations show that the function T is smooth while the 
function \1/ has a discontinuity at the zero time. These results take place for q < q cr , including 
q = 0. For q = one can see that an assumption about the parity properties implies that \1/ 
in fact has a discontinuity at t — 0. Indeed, equations of motion (j2.22j) can be rewritten as 
a system of equations on semi-axis: 



(3.26a) 
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(C_tf)(t) = T(t)tf(t), 0<t<oo, 
where C + and C_ have kernels if + and K_, respectively, 

K + (t,t') = K(t-t')+K(t + t'), 



(3.26b) 
(3.27a) 



Ji_ (i , t') = AT(t - if) - K{t + If) 



(3.27b) 



if(t) = -^exp(— t 2 ). These kernels are positive when both t, if > 0. If one assumes that 
< ty(t) < oo for < t < oo than it follows from the equation 



that T(t) > for < t < oo. This guaranties that (CT)(t) > for < t < oo. According 
to ()3.2fiaj) this means that * 2 > for t > 0, i.e. \1> 2 (0) > 0. Since \l> is an odd function it 
has to have discontinuity at zero time. 

This discontinuity disappears when we perform smoothing of this function with the kernel 
K to obtain the function ip. In fact this function represents the time evolution of the tachyon 
field. Shown in thin line, for comparison is ip a pp defined as smoothing of the solution to 
equation (|2~23|) . 

We can consider a function ip(t) starting at zero time from false vacuum and rolling to 
the true vacuum at large times. For q = this is a monotonic function. As it was shown in 
j2H] for Q 2 < Qcr there are small oscillations with decreasing amplitudes around the solutions 
to equation (|2.25jl . The presence of this oscillations (in fact oscillations with a complex 
frequency) results from the presence of infinite number of derivatives in a linear equation 
describing a deviation of ip{t) from ip q= o(t). 

4 Energy and Pressure 
4.1 Conservation of energy 

In this section we derive functional expressions for energy for tachyon equations (j2.18j) and 
for approximate equation ()2.23j) . Let us consider the system of integral equations (j2.19jl . The 
conserved energy is a sum of three terms representing kinetic energy E^, potential energy 
E„ and non-local terms 



T 



(3.28) 



E{t) = E k {t) + E p {t) + E nl {t), 



(4.29) 



where 



E k (t) = \{d^)\ 



(4.30) 




(4.31) 
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E nl (t) = lj o dp (e^T*) V (d t e E * 92 y) + 1 jT dp (e"^ 2 vl> 2 ) V (de^x) , (4.32) 

here and below d = f t and A d B = AdB - (dA)B. This expression is analogous to 
an expression for the energy obtained for p-adic equation from calculation of the energy- 
momentum tensor ^01 EE] • 

One can immediately see that this energy conserves. Indeed, 

= --vdv + q 2 (di))(d 2 i)) - i)di) + -(dT)V 2 + T^d^f 
dij 2 2 

+- f 1 dp(e- E92 T^(de Ed2 ^) + - I' dp(e- E92 y 2 )l?(e Ed2 dr). 
8 Jo 16 Jo 

Using the identity 



- m / dp(e- mpd2 (p) d 2 (e- m(1 - p)a V) = ^e" ma2 0, (4.33) 



we get 



= --vdv + g 2 (<9?/0(<9 2 ?/0 - W + -(<9T)^ 2 + T^<9^ 
2 2 



= ^8v[-v + e"^ 2 (^ 2 )] + d^[{q 2 d 2 - + e"^ 2 (T^)] = 0. 

On solutions to equation of motion the nonlocal term can be presented as the sum of the 
following terms 

E nl (t) = E nll (t) + E nl2 (t) + E nl3 {t) + E nl4 (t), (4.34) 



where 



E nll (t) = ~j\p(e^ d2 (- q *d 2 + l)ip) (e-^ 2 a>), (4.35a) 
E nl2 {t) = - l -j\p(e^ d2 {-q 2 d 2 + l)d^(e-^ d2 d^), (4.35b) 
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E nl3 (t) = dp (e^v) [e-^ d2 d 2 v) , (4.35c) 

E nU (t) = ~J dp(e^ d2 dv) (e-^dv). (4.35d) 

It is interesting to compare the expression (|4.34j) for the nonlocal part of energy with the 
nonlocal term corresponding to the energy for equation (J2.23|) . The latter has the form (we 
use calligraphic letters for quantities defined on solutions to equation (|2.23J1 ) 

S(t)=S p + £ k + S nl (t), (4.36) 
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where 



4 = fwf, E p = -~i> 2 + ^\ (4.37) 



and the nonlocal term mentioned above is given by 

E nl (t) = 1^ dp (e^ a V) *d (de***) . (4.38) 

Here ^ = ^ app and ip = ip app . This nonlocal term takes the following form on equation of 
motion 

£ni(t) = S nll (t) + = I J dp [e^ 9 \-q 2 d 2 + lty) V (ae-^ a V) . (4.39) 

Comparing (|4.39|) to (J4.35aj) and ()4.35b|) we see that functional expressions for S n n and £ n i 2 
are the same as the ones to E n u and E n i 2 . 

In Fig |U we plot functions £ p (t) + £k(t) and £ n i(t). We see that up to inaccuracy in 
numerical calculations the energy is conserved. To calculate the nonlinear terms numerically 
we actually use another representation for £ n i(t), 

E nl (t) = -J dp(e^ d \-q 2 d 2 + l)^ V(<9e^). (4.40) 

f!4.40|) allows us to use the integral form 

£m(t) = i jf dp{K x y){K 2 *) - i jT dp^)^*), (4.41) 

where #i = (-g 2 <9 2 + 1)C^, K 2 = <9 2 C|, K 3 = (-q 2 d 3 + d)C^, K A = OCe. In actual 
calculations we take e ~ 0.02. 




Figure 4: Testing the constancy of the energy for the interpolating solution to equation (|2.24|) . 
This figure shows £ p (t) + and £ n i(t), calculated from (|4.37j) and (|4.41|) . respectively. 
Their sum is equal to zero up to numerical errors. 
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4.2 Pressure 

The pressure is defined in terms of the energy-momentum tensor P(t)i = —T- (no summa- 
tion). Since we consider homogeneous configurations we will omit the vector index %. The 
pressure has the form 

Pit) = —E + 2E k (t) + 2E nl2 {t) + 2E nU (t), (4.42) 

or explicitly 



1 r 1 

Pit) = -E + q 2 (d^) 2 - - A I rfp(eV 92 (- g 2 9 2 + l)^)( e -V 92 ^) 

(4.43) 

dp(e~^ a dv)(e s Ed dv). 
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The analogous formula takes place for bosonic SFT [TBI I2H1- For equation (|2.24|) the 

pressure P is 

Pit) = -E + 2£ k {t) + 2£ nl2 {t), (4.44) 
i.e. the functional expression for this pressure is given by 



1 



V(t) = —E + q 2 {d^)) 2 - - / dp{e^r a \-q 2 d 2 + l)^)( e -T £< 3 d t V>)- (4.45) 

4 Jo 

For a large time the pressure in both cases is equal to —E. 

One can see that the pressure (|4.45jl at the zero time is equal to the total energy E, 
i.e. P(0) = E. This can be shown by taking into account the explicit relation between the 
energy and the pressure ()4.44|) and the fact that for our solution S n n(0) = 0. Indeed, 

V(0) = —E + q 2 (diP) 2 (0) + S nl2 (0) = —E + g 2 (c>^) 2 (0) + 2E- q 2 (diP) 2 (0) = E. 

For the case of system of equations (J2.19J) the pressure at t = is given by 

P(0) = -E + g 2 (cty>) 2 (0) --J dp (e^ d \-q 2 d 2 + 1)5^) (e^V 92 ^) . (4.46) 

One sees that this functional expression coincides with the expression (J4.45|) that gives the 
pressure for equation (|2.24|) . In FigEl we show the pressure Pit) for equation (|2.24|) as a 
function of time. To obtain this picture we use 

1 r 1 

<p( t ) = -E + q 2 (dij) 2 (0)--J dp(K 3 V)(K^). (4.47) 

The minimum of the pressure is at t=0. Note that for large times the value of P is positive 
and equal to — E = 1/4. 

In the above consideration we take into account only the energy of tachyon field. Follow- 
ing the Sen conjecture one has to add to the action an extra term representing the D-brane 
tension. This term shifts on a constant the tachyon potential making it equal to zero at its 
minimum. This constant also shift the pressure [2H] and we get that for our solutions the 
pressure goes to zero for large times. 
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b 



Figure 5: Calculating the pressure on the solutions FigG]to equation (|2.24jl : a) for q = 0, 
b) for q 2 = 0.96. This figure shows V(t) has minimum at t — and goes to 1/4 at the ±oo 

5 Conclusion 

In this paper we have studied the SSFT tachyon classical time dependent equations of motion 
describing non-stability of non-BPS D-branes. We have found the special time dependent 
solution to this equation that describes the rolling tachyon. It starts from the perturbative 
vacuum and approaches one of stable vacua in infinite time. We have argued that the 
pressure on this solution tends to zero for large times. These results have been obtained in 
the first non-trivial level within the level truncation method. In the next paper we will take 
into account higher level fields including vector fields and study a dynamics of a D-brane 
tension. 
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